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ABSTKACT 


vl/ 

Spatially  dependant  autoregressive  models  in  m dimensions  ar« 
defined.  The  oontitlona  for  etatlonarity  and  invertibility  are  de- 
termined, The  autocorrelation  function  and  Yule-Walker  equations 
are  obtained  for  the  general  case^  and  as  particular  cases 
z(t)  - f(x^,  x^;  t),  for  -oa  < t^,  for  special  discrete  values 

Xj^i»  and  t^,  for  various  grids  in  (x^,  x^)  plane  and  for  orders 

1 and  2 in  time.  The  spectra  are  obtained  for  these  pax*tlcular 
cases,  and  same  results  for  the  partial  autocorrelation  function. 

All  results  are  new.  The  notation,  definitions  and  assumptions  are 

those  given  by  Voss  et  al  (1977 )•  We  assume  stationarlty  of  z . 

x,t 

over  time  t,  where  x - (x^^,  x^,  ,,,,  x^)  an  m dimensional  vector. 

We  assume  the  covariance  structure  as  given  by  Hannan  (1970),  with 
C ^>0,  and  all  covariances  existing,  Non-statlonajy  models  will 
be  considered  in  later  papers. 


i 


W 


2 


1«  Dei'lnltions,  Model,  Assumptions,  Conditions. 

The  notation,  definitions  and  assuirptlons  are  those  given  by  Voss 

et  al  (1977).  We  assume  stationarity  of  z . over  time  t,  whore 

x,t 

X - (z^,  x^,  ...  , x^)  an  m dimensional  vector.  We  assume  the  covariance 

2 

structure  as  given  by  Hannan  (1970),  with  OT  >0,  and  all  covariances 

2 

existing.  Non-stationary  models  will  be  considered  in  later  papers. 

The  autoregressive  model  in  m dimensions  is  defined  as 

(1.1)  ^Xft  “^n-  -#»^k-l  i^n,k  *x+n,  t-k  * *x,t  * 

In  (1.1)  x^,  x^,  ...  , x^  are  m space  variables,  €uid  t is  the  time 

variable.  The  a's  are  independent  shocks,  i.e.,  independent  random 

variables  with  zero  mean  and  variance  ofia^.l^fiT'^ja  .is 

^ x,tv  a x,t 

independent  of  z . . except  when  k • 0.  The  autocorrelation  function 

X,  t“K 

of  the  whit,  noise  process  is  ,i«n  by  ' 0 WO 

for  - od<t<®*’.  We  assume  that  T . is  a weakly  stationary  process. 

The  autoregressive  process  can  be  thought  of  as  the  output  z . from  a 

x,t 

linear  filter  with  transfer  function  4*  ”^(B„,B. ),  iidian  the  input  is 

X w 

wliite  noise  a . In  this  paper  we  are  interested  in  some  special  cases 
x,t 

of  (1.1)  in  which  only  a finite  number  of  coefficients  are  non>zero,  i.e.. 


(1.2)  z^^^  ■ ^ ^ n,k  *x+n,  t-k  * *x,t  * 

where  p - (Pp  p,^,  ...  , p^^^),  and  q - (q^^,  q^,  ...  , q^).  In  this  case 

^(B  .BJ  is  finite.  Therefore  for  invertiblllty,  no  restrictions  are 

X Xr 

needed  on  the  parameters  The  autocovarlanoe  function  of  Uio  auto- 

regressive  process,  AR,  is  obtained  by  multiplying  throughout  by 

and  then  by  taking  expectations.  We  can  obtain  autocovariances  via  the 

oo  ^ oO  c d 

(B,B.)“iL 

- ~ X t C"- oo  (!*•-<»  cd  * 
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for 


where  is  the  coefficient  of  both  and  B”*^  B”^  , and  y^  is 

the  variance  of  the  process.  It  can  be  shown  that 

(1.3)  r - Crl  1"(Bj^,B^) 

where  V(B^»B^)  - 4’“^(B^,B^)  and  t(F^,P^)  - 4»"^(F^,P^). 

1.2 

Spectrum.  If  we  substitute  B.  e”  and  B 

1 - j - W ind  i ■ rTi  in  (1.3),  we  obtain  the  power  spectrum.  Thus 
spectrum  of  AR  process  is 

(l.h)  P(f,g)  - 2 <rl  1 4^(0'^^^^  , e"^^”""^)  I , 0 ^ f - 

0 - gj  - ^-S  , g - (g]^,  g2»  •••  » 8jj)  i 4*  (e“^^*,  - 

f>(e-^2irg,  ^ ^.i2irg^^  ^ ..i2rg^^  ^.i2trf 

Stationarity.  The  stationartly  conditions  for  AR  process  are  generaiiaations 
of  Box  and  Jenkins,  i.e.,  for  stationarity  4^  ”^(B^,B. ) must  converge  for 
some  sets  of  values  of  B ,B. . These  conditions  are  obtained  for  various 

X Xt 

models  in  the  next  two  sections. 

In  section  2 we  consider  the  case  of  m ■ 1,  and  In  section  3 we  con- 
sider the  case  of  m ■ 2.  In  each  of  these  sections,  various  models  are 
discussed  discussed  in  detail.  All  zero  dimensional  results  of  Box  and 
Jenkins  (1970)  are  special  cases  of  the  more  general  models  in  this  paper. 

2.  Autoregressive  Models  in  m - 1 dimensions. 

For  m ■ 1,  (1.2)  can  be  written  as  (writing  x for  x^  and  n for  n^) 

(2.1)  ^x,t  ” ^n“.-p,  ^ k-1  ^n,k  *x+n,t-k  * ®x,t. 

The  process  defined  by  (2.1)  is  denoted  by  AR  ( r;  p^^,  q^^)  , and  can  be 
wrttten  as, 

C2.2)  «|>  \ ^ , -here 

(2.3)  4>  (B,,B^)  . X - T„,k  K • 

The  autocorrelation  function,  Yule-Walker  equations  and  stationarity 
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conditions  are  obtained  for  various  special  cases  of  AR  (r;  p^,  q^)* 
Since  the  series  4’(b^.bj  ln(2.3)  Is  finite,  no  restrictions  are 

X w 

required  to  ensure  invertibility. 


2.1  The  Model  AR  (1;  1,0) 


writing  < 

Pi  I" 

\,t  • 'T’A.t-i  ' TsS-i.t-i  * Vt ' 

can  be  written  as  • *x,t' 

(2.1.2)  4>i(B^,Bj)  - 1 B^_. 

Autocorrelation  Function,  Multiplyinc  throughout  in  (2.1,1)  by 

(2.1.3)  *x_in,t-n®x,t  ” Tl^x-m, t-n*x, t-1  *^2*x-ni,t^  *x-l,t-l  * *x-^B,t-n  *x,t  * 

Qa  taking  expected  values  In  (2.1.3),  we  obtain 

ya,n  " f lY«,n.l  " T2Y„^l,n-l 

m - 0,  n ? 1;  m>l,  n ■ 0;  and  m - 1,  n ^ 1. 

(2.1.5)  Also  Yoo  ■ firol  * ^2  Yu  * Yo-1  *Yol  ““  Y-1-1  'Yu- 

Ol  divldijig  by  in  (B.l.B)  Mid  (2.1.5). 

/°»,n  ■ f 1 f m,„.l  *f2  ^».l,iv.l> 

not  m - 1,  n = 0;  and 

(2.1.7) 

From  (2.1.U)  and  (2.1.o),  we  see  that  both  autocorrelation  and  auto- 


covariance  functions  satisfy  the  sane  form  of  difference  equations,  Ftom 


(2,1,7)  , the  variance  (T  may  bo  widtten  as, 

2 


(2.1,0)  (Tl . 0-1/  (1  -<p^  />u)  . 

Prom  (2.1.2)4>^(B^,B^)  - 1 -(<}>i*f2®x)  ®t  «ndS'(B^,B^)  -4*^ 


1 
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I 

I 

I 
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} 


If  we  suDstitute 


-2Trif  . ^ 
e and 


e”^  ^ , i "TT^  in  tlie  auto- 


't  ■ X 

covariance  generating  function  (2*1«8),  we  obtain  the  power 

apectrom.  Thus  the  spectrum  of  AR  (I5  1,0)  is  given  by 

(2.i.9)  p(f,8)  ■ ■ 

icrl/  +^2  ■*■^^2  - 2^<p^ooe  2irf  ♦'^2®°®  2ir(f+g^J  , 

0-g-^ 

Again  when  (^2  " ^ Sot  tlio  corresponding  result  for  zero  ditnanslons. 

Statioiiarity.  For  stationarity,  the  autocovariances  and  autocorrelations 
must  satisfy  a set  of  conditions  to  ensure  stationarity*  The  conditions 
can  be  combined  in  a single  condition  that  the  generating  function  ^(B  ,B. ) 
must  converge  for  )b^)  - 1 and  | B^l  - 1*  See  Theorem  (2.1,12)  below.  Now 

C2.1.10)  . (1  - bJ-1  ('Ti*T2V‘‘  ^ * 

From  (2.1*10)  we  see  that  for  stationarity  ^ 2®x ^ ^ ^ 

Taking  - 1 and  B^  - -1,  the  parameters  of  AR  (1|  1,0)  must  satisfy 
1^1*  <^2^"^^*  and  | -(p2l<l.  to  ensure  statlonazdty.  The  two 

conditions  can  be  combined  into  a single  condition* 

(2.1.11)  |<y^l  ♦1^2!  ^ 

This  region  is  shown  in  Fig*  1 


■fr 


We  briefly  give  a justification  for  Ib^  I - 1.  Operationally  tlie  condition 
for  stationarity  may  bo  written  as 

■ * ‘fl  *f2“x>  “t  ‘ (fl  bJ  ^ 
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“ ‘x,!  " <Tl^f2®x>  S,t-1  " ^,1^2  " — 

" •x,t*'^'l*x»t-l*?2*x-l,t-l  *fl*x,t-2*  ^'^l‘^2*x-l,fc-2+f2a^_2,t_2^  ••• 
This  series  does  not  involve  B. • Hence  it  is  zero  dimensional  in  x. 

So,  |b^)  - 1,  as  must  be  tfio  case  for  stationarlty  for  any  B operator 
since  we  assume  stationarity  in  -^x,  tj^  , x a vector* 

(2.1*12)  Generalization  of  the  theorem  in  zero  dimensions* 

The  roots  of  tiie  characterif*-ic  equation  must  lie  outside  the 
unit  circle  in  each  variable,  B^  when  all  other  B's  ax^  set  to  one* 

The  condition  must  also  hold,  as  stationarlty  in  n dimensions,  m 4-  1 
variables,  requires  stationarlty  in  every  direction  in  the  m 1 
variables*  The  same  condition  holds  in  case  of  Invortlblllty  for  the 
Ml  models. 

Autoregressive  paramo ters  in  terns  of  autocorrelations  Ikile-Wivlker  oquatlons* 

From  (2*1*6),  we  obtain 

(2*1*13)  /’oi  - fl  "?2  riO 

Pu.  " f 1 /^lO  * ? 2 Hence 

(2.1.1U)  <^1  ■ (/oi  " Fu.Py^ / ^ "P y? 

‘f  2 “ “ foi  Py?  / ^ ^ “P 10^ 

Again  from  (2*1*6)  we  obtain 

(2.1.15)  • f 1 ^o,k-l  * T2  fl.kJ  > 

(2.1.16  /’,,o-flfk,l  ‘fkPk-l,!  • 

(2.1.17)  - cj)j^  /’k,k-l  * f 2 /’k-l,k.l’ 

PonnilaB  (2.1.7),  (2.1.13),  (2,1.11|),  (2.1.15)  and  (2.1.16)  can 
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be  used  recursively  to  compute  various  autocorrelations 
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A special  c/isu  of  occutg  if  “ /^U  * (i^.l.lU) 

‘92  ‘ ^01^  ^ Pio  ^ As  an  ox.'unpl6,  if 

/^oi  Vn  "^10  • 9l  ""^2  “ 

Another  special  case  occurs  if  ■ 1,  then  from  (2.1.11)  we  have 
^ ♦<^2*  /^ii  ■ 9i  ^ f 2 meaningless  unless  “/^11» 

and  of  course  “ /^Ol  " /^ll  * means  '*’92^  ^ 

one  ccnditlon  as  alreac^  seen  in  (2«l»l0)  for  inTortibillly*  Such  cases 

2 

hoirever  are  singular  since  we  insist  ^ ]^q  ^ 

Further  wc  note  that  <p  and  <p2  detenninod  by  /'lO*  C YL' 

Conversely  if  and  arc  chosen  in  accordance  with  (2*1*13))  then  given 

any  one  of  q-^*  fvj*  fll*  other  two  are  det<jnninod, 

(2.1.18)  As  an  example  wo  take  cp^  ■■  •2>  cp^  " » 

‘xt  ■ *2  »X,V1  - 'x.l,t-l  * ‘xt  • 

The  sot  of  values  for  the  ^ 's  is  determined  by  the  equations  (2.1.3)> 
given  <?p2  l^iQ  ^ solve  for  and  ^ jj*  We  choose  - .5 

and  consequently  » /^ll  “ “•^*  autocorrelation  function 

is  given  by  (2.1.6),  (2.1.15),  (2.1.16),  and  (2.1.17).  In  all  cases 

“"/^no*  the  special  case  ■ .2,  Cp2  “ ti«  auto- 
correlation function  jO is  given  by  the  diagram  whero  tlie  x axis 

f mn 

roprosents  m and  the  y axis  n. 
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m 


.5 

.25 

.125 

.0625 

1 

-.5 

-.25 

-.125 

-.0625 

28 

-.Oii 

.25 

.125 

.0625 

08 

-.18 

.07U 

-.125 

131 

-.083 

.120U 

076 

-.095 

072 

.03125 

-.03125 


n 
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The  analysis  for  the  model  AR  (1;  0,1)  Is  given  by 

(2.1.19)  Irt  ■ fi\,ul  . tsVl.t-l  * ‘xt* 

It  is  similar  to  model  AR(l;  1,0) » The  results  are  analogous  if 

we  replace  the  parameter  2 Tl  operator  byP^. 

2.2  Tlw  Hfxtel  AA  (2}  1.0).  AR  (2;1.0)  Is  glvon  by 

(2.2.1)  *3jt  * Tol  ®x,t-l  ^ '^llVl,t-l  * T02*x,<i-2  ^ 1^12®x-l,t-2  * *xt  * 
(2-2.2)  ‘♦’2  (VV*xt'‘xt 

(2.2.3)  4>2  (B,,B,)  ■ 1 - ( <?01  ‘Tu'x  * ^02^  ’ 'fl2Vt)  “f 

Conditions  satisfied  by  tiio  autocorrelations.  The  correlation  matrloes 
for  the  Yule-Walker  equations  in  t and  x are, 

(2.2.1*) 


or 


0*2  p 


t X 


<r‘ 


1 r 

rio 


10 

1 


and  Cr^ 

z t 


^ Tul  fo2 
/’oi  ^ /’oij 
L/^02  fbl  ^ 


p 

Prom  (2.2.I4),  positive  definiteness  implies  1 “/^jj  ^ 0 , 


" hi  fn 

poz 

fl2 

^ fol  f 02 

hi  ^ Ao 

Pol 

fn 

(2.2.5) 

foil  ^ C 01 

> 0,  and 

f\l  f 10  ^ 

Cn 

P 01 

poz  foi  ^ 

f^02  Pol  P 11 

1 

Pio 

Pl2  Pu  Pol 

P LO 

1 

Autocorrelation  function. 

Multiplying  (2.2.1)  by  z . , we  got 

x-n,  b**n 

sv  ^ 

*x-m,t-n*xt 

Toi*x-m,  t-n®x,  t-1 

*Pii*x^^t_n®x-l,t-l  * ^02*: 

z 

X,  U2  x-m,t-n 

>0, 


* Tl2®x-m,t-n*x-l,t-2  * ®x-m,t-n  ®xt. 

Taking  ex.ected  values  on  both  sides,  we  gut 

‘2-2-M  ■ <^01  Y.,.^1  * ‘fuYs^l.^-l  * toe  Ym.n-2  * Tl2  Ym.1,.^2 

for  m ■ O,  n - 1;  m > 1,  n ■ 0;  n - 1,  m ^ 1. 


Also 
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roo''Poirol*1'lliril*<?02r02*'Tli^l2* 

Dividing  by  in  (2*2«6)  and  (2.2a7) 

/>„  • tol/"n,n.l  * Tu /’.-l.n-l  * foyA.n-y  * Tl2  /’»-l,n-2 
for  m - 0,  n ^ 1;  m>l,  n - 0;  ^ ^ 1,  n - 1,  And 

(2.2.9)  1 - ‘ Tu/’U  • To2/>02  ‘ tl2/’l2  * '^lAl  ' 


The  variance  from  (2.2.9)  is  given  by. 


(2.2.10)0'^  - "Tol/^01  ■ “ f 02/^02  “ ‘fl2/’l2^. 

Spoctnao.  Let  H^(B  ,B.)  ,B.)  , and  from  (2,2.3) 

X V*  c X w 


t j ■ I - ( 1>ol  • f ll“x  * ‘t'02>t  * 'Pl2Vt)  ®t  • 

.-2£irf  „ . -2iire 


Substituting  B. 


and  3 ■ e 

X 


in  the  autocovariance 


generating  function  f^(B^>B. ) given  in  (1.3) » we  get  the  power 

X w 

spectrum.  Thus  the  spectrum  of  AR  (2;  1,0)  is  given  by 


(2.2.11)  , 

P(f.8)  ■ 2U^/  (i  -foi'*'*'” 


11' 


,-2OT  -2irg 

Stationarity.  As  in  section  2.1,  for  stationarity  '4^(B  ,B.)"  *^”^(B^,B  ) 

X V Xu 

• - <foi  ' foy®!  * '?12Vt>  ]■" 

■ ^d-0  'fol  * Toy  ^ Tu®!  * TlyBgBt  ) 3 °t 

must  converge  for  | - 1 and  |b^|  ^ 1 . 

This  implies  that  the  roots  of  1 - ((<^q-^^  *^ll®x^^t  * ^To2  ^t.)  " ^ 

must  lie  outside  the  region  (b^|  - 1 and  |b^|  'i'  1,  The  conditions  are 

(2.2.12)  ^02*  '?12*?01^fll 
t02*  Ti2"1’01-  TU 
02  “ ^^12  Tol  ■ ‘Til  ^ ^ 


To2'?12  -?01"Ti1 

Ifoyl  MTiyl  <'• 


< 1 


and 


M 
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Yule-Vfalker  Equations. 

(fcltiplyii*  li.  by  Vl,t-2 

respectively  and  taking  expected  values  on  both  sides,  we  get, 

(i:.2.1ii)  ♦ fllPlO  ^ 1^02  * ‘^12/^11  * 

(2.2.15)  /^ii  ■ 9oi /^lO  ^ 'f  11  * ?02 /°11  * ^12^01  * 

(2.2.16)  “ 'f'oi  ^01  '^11  <°11  * ^*02  Ti2  PlO  * 

(2.2.17)  /^i2  " ‘^01 /^ll  * ^11^01  * '?02/°10  * 

These  equations  can  be  solved  and  parameters  '^^y*  ^02*  ^12 

can  be  expressed  in  terms  of  autocorrelations.  Also  from  (2. 2.6)  we  obtain 

(2.2.18)  p Ok  ■ ?01  D,k-1  * ^uf l,k-l  * To2  /^0,k-2  * f 12  Pl,k~2  * 

(2.2.19)  " 'f  Qi  /^k-1,1  * ^ 02  /^k2  ■*‘1^12  fk-1,1  * 

(2.2.20)  • foif^k,k.l  * fu/^ k-l,k-l  ^ 9o2Pk,k-2  * ^ 12  /’k-;L,k-2  » 
where  k > 1.  Formulas  (2.2.18),  (2.2.19),  (2.2.20),  and  (2.2.8) 

can  be  used  recursively  to  compute  various  autocorrelations. 

2.3  Partial  autocorrelation  function  for  AR  (1;  1,0).  The  partial 

autocorrelation  is  a device  which  helps  us  decide  which  ordor  autoregressive 

process  to  fit.  irfith  the  models  given  in  sections  2.1  and  2.2,  we  prove 

the  following  rcsvilt. 

Theorem.  If  the  tiue  model  is  AR  (1;  1,0)  , then  - 0 and  ^ » 

assuming  ft 

Proof.  From  (2.2.1U)  and  (2.2.15)  , 

/^oi  "/^ii  " ^oi^/°oi  ~ f^ioFu.^  * lofol  ~Pu)  *1’o2^  foi  “/^ll  ^ 

Substituting  tiie  values  of  and  ^11  “ ^2 

(2.1.12)  and  after  simplification,  we  get  ^q2( /^q1  “ 1 ~Piq)  - 0. 


But  1 - / by  (i;,<;.bi)  and  /O  by  aBOumptlon,  tharuloro 

■ 0.  Similarly  it  c:in  be  shown  that  ■ 0,  which  compietefl  tlw 

proof#  It  should  bo  noticed  ttiat  all  other  parameters  like 
and  for  k - 2 


are  zero* 


2.1  The  Model  AH  (1;  1,1) 


'xt  " *^01*x,t-l  ^ ^-ll\-l,t-l  * '^*ll*X-*'l,t-l  ' “xt 


■*•  a . , which  reduces 


Hic  model  AH(1;  1,1)  is  (;iven  by 
L*: 

inmediatoiy  to 

* V =!,«,  -T^ 

“ /°-u  ■ ^11  • 

(2.a.2)  Hence  ♦ ■ 1 - B^. 

The  autocorrelation  fixnction  is  found  by  multiplying  (2.U#1)  by 
to  got. 


r' 

z 


x-m,t-n 

^x^.t-n\t“  '^OlVm,t-n^,t-l 

* z n.  a . We  take  oxpocted  values  to  got 

x-m,t-n  xt# 

(2,h»h)  Y ,n,n  " ^OlYm,n-l  * ^ll^Vm-l,n-l  *ym+l,n-l^  * 

(2.I..5)  y„,  ■ <Pol  roa  * filYu  * '^a  > 

Pn,„  ■ ‘Pol  ^,.^1  * 

^ 't 

m>l,  n -0;  m - 1,  n - 1#  From  (2.1uf>)  if  wo  divide  by  (T^* 

(2.1,.7) 

-I’lIT  f 

The  power  spoctnom  ifi  found  as  usual  by  substitution  of  ■ e , 


B 


, . ,2lTg 

’ X ’ 


.f:r 


in  the  autocovariance  gonoralin*' 


function  r(B^,P^,B^) 

(2,h.d)  p(f,g)  - 2<r^  |l  - <p^(  , 

0if5ii,0ig^^# 

The  conditioiui  for  stationnrity  an,  found  by  Bul.tlnu  U - <1,  F - ♦ 1, 
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and  B.*  - 1 using  the  same  arguments  as  before.  Now  '^'(B  ,F  -B  ,F  ,B.  ) 

V X X w X X V 

■ - (Toi*Tu<V^x)  ■ K (Toi“fu  ‘V*'.>  )“  4 • 

Hence  tlie  restrictions  on  the  parameters  aret  Itol* 

(2,k»9)  or  ^ ^ tiie  results  given  in  (2,1), 

Tiie  Yule-Walker  equations  are  found  from  (2,/<,6)  using  m»0,  n*l,  and 
m = n ••  li 

(2.U.10)  /'ll  "fol/'lo  fll<^Y2U^*  Consequently 

■ (^02  ^^■*’/’  20^  " ^/’lo/ll)  / (}  * f 20  ' ^f^io) 

?11  “ C/'ll  ‘ f 01  Pio  ) / ( ^ ^/'20  "/'lo}  • 

The  x-ocursion  formulas  for  the  autocorrelation  fuiiction  axt:  » 

(2.U.12)  ^ '■Toi/'o,k-l  * ^Tll/’ l,k-l  ' /icO  "TcJl/’ k,l  ^Tll^/’ k-l.l^/’k+l,!^  * 
/^kk  “fol/'k^k-l  * ^ll^/’k-l,k-l  /°k+l,k-l^  , k i 1. 

Now  - 0,  if  -./’oi  Ao*  ^^tii  " ‘ foiCio 

O 

If  however  1 + " 2 ■ 0,  tlion  as  a consequence  ■ f^olf^lo  * 

and  then  - 0,  The  condition  1 ♦ “ 2 'ir  0 is  satisflod  wtion 

m - 0,  or  “ 0. 

3*  Autoregressive  Models  in  M ••  2 Dimoxisions 
For  m ■ 2,  (1«1)  can  be  wid.tton  as 
fin'**  ■T'**'  ^ 

' * ' *x^,x^,t  ^ n^“- t»»  ^n2*-*“^k-l ‘fnj,n2»k*Xj^+nj^,X2'*’n,^, t-k  * ®x^,x,^,t  ' 

Consider  a special  case  of  (3»1)  in  which  only  a finite  number  of t^'s 


are  non-zero,  i,e.,  m-  2,  in  (1.2)  , 


(3.2) 


\,X2,t-^x^-p^In^-P2^Jnl^n^,n2,kV^’V'^»^‘‘  * « 

The  autoregressive  process  defined  by  (3.2)  is  denoted  as  Al!(r;Pj^,qj^;P2»q2^  » 

and  can  be  written  as  T,(B  ,B  ,B. ) 'z  . ■ a . 

3'  x^*  x^*  t^  x^,X2,t  x^,x^,t 


» 


whore 
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0.3)  4>,(b^,b,^,b,)  ■ 1 

Autocorrelation  function,  stationarity  conditions  and  Yule-Wd.  ker 
equations  are  obtained  for  various  special  cases  of  AR(  r; 

Clearly  various  results  obtained  reduce  to  the  corrospondin^j  results  for 
m •>  0 or  1. 

3.1  The  Model  AH  (1;  1,0;  1,0).  

Prom  (3.2),  AH  (1;  1,0;  1,0)  is 

*xyt  “^001  *xy,t-l  Tl01*x-l,y,t-l  * f Oll*x,y-l,t-l 
• f * “xyt  > by  X and  b,  y for 

conyenienoe  and  alao  f.jpiby  f o.u  Toll* 

First  we  consider  two  interesting  special  cases  of  (3*1.1) 


Special  Case  1. 


■ 0,  i*e«,  we  consider  the  model 


(3.1.2)  x^  “Tool*xy,t-l  * Tlol®x-l,y.t-l  * T 011*x,y-l,t-l  * “xyt 

(3.1.3)  <t>3i(B^,Bj^,Bp  ■ a^  , wham 

(J.l.h)  4>  3i(B,,B^,Bp  . 1 - epppj  a <f  B^)  B^  . 

Autocorrelation  function.  Multiplying  in  (3,1,2)  by  z 0 ,r  « + r. 


taking;  the  expected  values,  wo  obtain 

(3.1.5)  y^nn  - ^ OOl'yLmyn-l  101  y)',-l,m,n-l  * ^OllTt  ,m-l,n-l  and 

at  least  one  of  w{.,m,n  greater  than  zero  and  others  zero  or  positive, 

Yooo  ■ TocaYool  * TioiYica.  “f  ouYou  * * 

2 

Prom  (3.1.5)  on  dividing  byyQ^  “(T^  , we  got  tlio  autocorrolution 
function 

(3.1.7)  PIj^  - * T 101  A -l,m,n-l  ^ Toil /’p,r..l,n-l  , at 

least  one  of  , m,  n ^0,  ni\d  otliors  zero  or  pnaltlvc.  Prom  (3.1.0), 
the  varifmcc  CT^  is  obtained  on  dividing  by 

(3.1.8)  (T^  ■ ■ f OOlYoOl  “ flOl  YlOl  " TollYol})* 
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Spectrum.  Let  ,B  ,B  ) - 4^"^(B  ,B  ,B  ) , and  from  (j.l.U), 

A y w»  X y w 

1 - < fool  ‘ floA  * foiA>  “t  ■ 


B «■  e 

X 


y'  t ■ I uui  I xui  X t u±i  y t t 

-2iTrg  ^ ^ autocovariance  generating  function. 


r-  f(B,.B  ,B^)  4'(F^  ,F  the  apoctrum  of  autoregreasive  model 


(3.1#2)  is  given  by 


0 = f^^-s,  V - ^ ^ , O^hf^s. 

Stationarity.  From  stationarity  converge  for  | j ^ 1, 

I B^l  t 1,  and  I B^l  ^ 1.  Prom  (3.1.ii) 

O.1.10)  t3i(B,.By,B,)  • (l  - .■p^^By)  B^)'' 

oa  , . 

^ d-0  ^fool  ‘*'*^lOl®x  ■^^011^  ®t  • ^ »10),  for  stationarity 

«-l-ll)l<Pool*  floA'Touh^  for|B^l  il,  (BjI  i 1.  Taking 

B ■ 1,  -1  and  B ■ 1,  -1,  thn  stationarity  conditions  nrei 
X y 

(3.1.12)  l^fooi  * Tioi  * Ton  I <1 

I Tool  “ flOl  Toll  I ^ ^ 

] T 001  * Tioi  "Ton  1 ^ ^ 

) Tool  " Tioi  "Toni  ^ • 

Yule-Walker  Equations.  Multiplying  throughout  by**  * 
and  y ill  (3.1.2)  and  taking  the  expected  values,  we  get 

(3.1.13)  /’ool  " Tool  Tioi  (^100  Ton  iio 
/’loi  “ Tool  /loo  ^TiQi  ♦ Ton^no 

/*^on  " Tool  /l)io  * Tioi  /iio  * Ton  . Let 

(3.i.ili)  1 f^ioo  no 

^ “ fioo  ^ Piio  - X p ^ " ^in  P P 

/ OlO  f li)0  / 110  M)lO'l(K)  -^ 

f JIO  f' LlO  ^ 


15 


(3.1.15) 

C 001 

/loo 

/olO 

LOl 

1 

/ no 

# 

/on 

/no 

1 

(3.1.16) 

1 

001 

010 

^ 101  ■ 

D"^ 

/l(X) 

/ lul 

/u.0 

^ 310 

/on 

1 

(3.1.17) 

1 

/ 100 

/1x)l 

Ton* 

D"^ 

/lOO 

1 

/ LOl 

/olO 

/no 

/on 

Special  Case  2 

• fin 

■ 0 end 

Tool 

series  in  m - 

L dimensions  and 

the  model  is 

0.  This  reduces  to  a time 


(3.1.18)  - ^l®x-l,y,t-l  "■^2*x,y-l,t-l  ^ ®xyt  * TlOl"  K*  ?0-ll  "^2 

(1 . - b.^y^) 

Autocorrelation  function*  Multiplying  throughout  in  (3*1.18)  by 


*x  £,y-ni  t-n  taking  expected  values,  wo  obtain 

(3.1.19)  ^ ^2ru^l,r^l  * Of4in,nr0,  and 

Yooo  • ^1  Yioi  * ^2  Yon  ‘’y  Togo  O.U20) 

(3.1.21)  ^ 1 * ^2 /i,n>-l,n-l  , and 

(3.1.22)  1 - /"lOl  ^2  /^on  * z * From  (3*1*22),  ttio 

variance  (T^  is  given  by  <J^  - Or^y(l  - - ^g/fn)  ' 

Spectrum*  From  (3*1*18) 

(3.1.23)  4>  (B^,By,B^)  - 1 - V 

(3.I.2I*)  Let  4'(B  ,B  ,B. ) ■ 4^“^(B  ,B  ,B. ).  The  autocovariance  generating 
3C  y x»  X y V 

function  is  given  by  T - 0'^'P(B^,By,,B^)'|'(F^,Fy,F^).  Ut  B^  - 

B ■ and  B ■ The  spectrum  of  the  AR  process  (3*l*ld) 

X y 

is  given  hy 
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(3.1.25)  ■ 2 T^^yji . 1 ' i 

oifi^a,  o<g£Jj,  o^hSh* 

Statiunarlty*  Expanding  (3*1«18)  as  In  (2.1.10),  the  stai.ionarlty  conditlone  are 
I ^l®x  ■*■  ^ ®xl  " and  I I ^ 1,  which  glvoa 

0,1»26)  I ^2  I ^ * I " ^2^  ^ conditions  aie  sijnilar  to  (2,1,10), 

TUle-Walker  Equations,  Multiplying  by?  - . . and?  and  taking 

x-i,y,t-i  x,y-i,t-a 


2 

expected  values  on  both  sides,  dividing  by  Yqqq  " 0"^  » we  get 

(3.1.27)  »d 

(3.1.28)  * ^2  * Solving,  we  get 

(3.1.29)  ^ j2J./llo)/t^  UO^ 

^ 2 “f/^on  ■ fioi  f no)/^  “/’llO 


This  model  and  all  results  obtained  are  similar  to  the  n ••  1 dimensional 
model  considered  in  section  2,1 


3.2  Other  Models.  Yule-Walker  equations,  autocorrelation  functions  and 
scxne  other  results  are  obtained  for  other  models  m ■ 2, but  are  not  reported 
in  this  paper.  Some  of  these  models  are: 

A.  AH  (1;  1,1, 1,1)  Nine  point  model  (10  parameters) 

B.  Special  case  of  A ^ point  model  (6  paz*ameter8) 

C.  With  B - AH  (2;1,1;  1,1) 

5 point  model  - (11  parameters)  • • 


0,  Special  case  of  AH  (1;,2,2;  2,2) 
13  point  model  (lit  parameters) 


• « 


B,  For  model  given  in  (3.1.1)  * 

U.  Conclusions.  The  general  autoragressive  models  in  m dimensions  have  been 
defined  and  general  theorems  obtained.  The  well-known  sero  dimensional  thooiy 
is  a sub-case.  InqDortant  cases  for  m ••  1 aid  2 liave  been  ccjnsidiired  and  tliuLr 
properties  investigated. 


' * 
r 
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